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1. Introduction 

Relativistic fluid dynamics has been applied with success to describe the 
dynamics of hot and dense matter created in relativistic heavy-ion collisions 
at RHIC, and recently, at LHC.^ Despite this success, the derivation of 
relativistic fluid dynamics from an underlying microscopic theory is still 
an open problem and is under intense theoretical investigation. This is 
mainly due to nontrivial consequences imposed by causality and stability 
in relativistic fluids.'^''' 
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In this paper we discuss a general method to compute all linear transport 
coefficients (such as the shear viscosity and its corresponding relaxation 
time) associated with the dissipative behavior of fluids. In particular, we 
prove that this prescription gives a value for the shear relaxation time that, 
under the 14-moment approximation, coincides with the one derived from 
kinetic theory. This shows that transient dynamics is determined by the 
slowest microscopic and not by the fastest fluid-dynamical time scale. Our 
approach^ depends only on generic analytical properties of retarded Green's 
functions. 



2. Definitions and Power-Counting Scheme 

Let us consider the following linear relation between a dissipative current 
J {X) and its thermodynamical force F {X) in Fourier space, 

J{Q) = Gr{Q)F{Q), (1) 

where we define the Fourier transformation in the following way, 

A{Q) = j <fXe-yci>{iQ-X)A{X) , (2) 

A{X) = j ^ cxp (-zQ . X) A (Q) . (3) 

Here, we adopt the following notation for the four-momentum and space- 
time coordinate, Q = {lj. q) and X = {t, x), respectively. In this article, we 
discuss particle diffusion, in which J is the diffusion current of a density 
n and F ^ dn, and shear flow, in which J is the shear stress tensor tt'^'^ 
and F the shear tensor a'^" = V^'^u'^^ with the fluid four- velocity u'^, 
y/i = Af'^d^, A^"" = gi^" - uV, and A'^^"''> = i(A'*«A''^ + A'^'^A''" - 
2A^-^A"0)A„0. 

Equation (1) contains all the information about the linearized underly- 
ing microscopic theory in the form of the retarded Green's function Gr{Q). 
When the system exhibits a clear separation between the typical micro- 
scopic and macroscopic scales, A and I. respectively, it is generally assumed 
that the dynamics of the system can be described in terms of a finite number 
of macroscopic variables. The macroscopic scale is associated with the varia- 
tions of the thermodynamic force F ~ and its derivatives d^F ~ 
The microscopic scale A is contained in the poles and derivatives of Gu{Q) 
and would be, for example, the mean-free path in dilute gases. The chal- 
lenge is to determine the equation of motion satisfied by the current J 
in terms of the microscopic information contained in the retarded Green's 
function Gr{Q) when such a separation of scales exists. 
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3. Gradient Expeinsion 

The separation between microscopic and macroscopic scales is normally 
used to determine the macroscopic dynamics of the dissipative current J in 
terms of an expansion in powers of Knudsen number, Kn = A/£ -C 1. The 
simplest realization of this idea is known as the gradient expansion.^ In the 
gradient expansion the dissipative currents are expressed solely in terms of 
gradients of fluid-dynamical variables. In practice, the gradient expansion 
is implemented by expanding the retarded Green's function, Gii{Q), in a 
Taylor series around the origin, Q = 0, 

Gr{Q) ~ Gr{0, 0) + d^GRiO, 0)uj + ..., (4) 

where the terms resulting from the expansion in q are not explicitly de- 
noted. This expansion leads to the following equation of motion for J 

J{X) = CF{X) + CidtF{X) + ... , (5) 
where C = Gr {0,0), Ci = du, Gr{Q) , .... We remark that the 

w.q— 

gradient expansion usually docs not involve time derivatives of the ther- 
modynamic force. This is not a problem since one can in principle use the 
conservation laws to replace time derivatives by spatial gradients. 

Apart from overall factors which ensure the correct dimensionality, C ~ 
A and Ci ^ A^, respectively. Thus, the term CF is of first order in Kn while 
C\dtF {t,x) is of second order in the Knudsen number, Kn^. If Kn ^ 1, 
higher-order gradients can in principle be omitted and one obtains a closed 
macroscopic theory for J. Navier-Stokes theory^ corresponds to tnmcating 
the series at first order while keeping higher-order corrections leads to the 
Burnett equations.'' 

This method, however, leads to unstable fluid-dynamical theories. In the 
non-relativistic case, the second-order truncation, i.e., the Burnett equa- 
tion, is unstable.® In the relativistic case, Navier-Stokes theory itself is 
unstable. ^'^ Recently, it has been shown that this instability in relativistic 
theories is intimately related to the violation of causality.^ Such problems 
prohibit the application of these theories to describe any realistic system. 



4. Computing the Relaxation Time from the Poles of Gr 

In the derivation of Eq. (5), two main assumptions were employed: the 
separation between the microscopic and macroscopic scales and the gradient 
expansion itself. The first assumption is essential for the coarse-graining 
procedure that defines fluid dynamics and cannot be disregarded. On the 
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other hand, the second assumption is solely based on the idea that the 
Taylor expansion is able to describe the retarded Green's function in the 
whole macroscopic/mesoscopic frequency domain. This need not be true, 
since it is known that Green's functions are not necessarily analytic in the 
whole complex plane. In fact, these functions can have singiilaritics in the 
complex w plane, which necessarily restricts the domain of validity of the 
Taylor expansion. 

In order to illustrate the physical meaning of such singularities, it is 
convenient to consider the generalized diffusion equation for J, also known 
as the Maxwell-Cattaneo equation.^ In this case, as in the Israel-Stewart 
theories of fluid dynamics, ^° the dissipative current satisfles a relaxation 
equation. In the linear domain, this can be written as, 

TR^J+J = DF. (6) 

In these theories there are two transport coeflacients: the relaxation time 
tr and the diffusion coefficient D. In Fourier space, Eq. (6) can be written 
in the form = Gr{uj)F{uj), where 

Gr{lo) = — \ . (7) 

TR W + i/tr 

It is important to notice the existence of a pole of Gr{u)) on the imaginary 
axis. The existence of this pole is directly related to the existence of a 
non-zero relaxation time and, consequently, to the causal structure of the 
theory.^ It is also clear from Eq. (6) that Pick's law and, consequently, the 
diffusion equation are recovered in the limit of vanishing relaxation time, 
tr ^ 0. Mathematically, this limit pushes the pole to infinity, rendering 
Gr{u}) an analytic function in the whole complex w plane. Only in this case 
a Taylor expansion near the origin has an infinite radius of convergence and 
is able to provide a good description of Gr{uj). 

However, in transient theories of fluid dynamics, where the relaxation 
time is necessarily nonzero, can a Taylor expansion still be used? Strictly 
speaking, the range of the Taylor expansion is limited by the existence of the 
pole. As was previously shown, the pole is directly related to the relaxation 
time and, consequently, the Taylor expansion is unable to probe time scales 
of the order of the relaxation time. 

5. The Role of the Analytical Structure of Gr ('*') 

Next, we consider a general analytical structure for Gr{Q). We assume that 
Gr{Q) has poles in the complex uj plane, wi(q), . . . , u;Ar(q), while the 
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structure in the q plane is analytic. Note that, in principle, the number of 
singularities does not need to be finite. Also, in many cases GniQ) contains 
branch cuts, especially in the UV region, oj/loi 1, which is sensitive to 
temperature-independent properties of the system. We only consider the 
region in the complex uj plane without branch cuts. Then, GniQ) can be 
expressed in the following way 

gr{q) = — "^'^^1 s = -f — / Ni ^'^r'^ ' 

^w-Wi(q) [w - wi(q)] • • • [w - WAr(q)J 

where fi{Q) and H!((5) are analytic functions in the complex w plane. We 
write the right-hand side of Eq. (8) in a polynomial form, 

l+$r(q) (-^^) + . . .+*.(q) {-^^f = {-if - ^i^j ' - - ^^(q)] , 

(9) 

where the Wi(q) are the zeros of the polynomial on the left-hand side and 
we introduced the functions $TO(q) 

$„(q) = Hr E ^. ' (a) " ^'^^ 

i<n.^„<jv'^^iW---'^«„.(q) 

Combining Eqs. (8) and (9), we obtain 

1 + *i(q) i-iLo) + ... + $^(q) i-iLof] Gr{Q) = ^"^^"^^^^.^ . . 

(11) 

After taking the inverse Fourier transform and expanding the functions 
^m(q) in a Taylor series around q = 0, we obtain a differential equation 
satisfied by Gii{X — Y). It should be clear that this will be a linear dif- 
ferential equation of order N in time. The equation of motion for J can 
be obtained by multiplying Eq. (11) by J and taking the inverse Fourier 
transform. See Ref.^ for details. The result is 

XNd^J + ... + xidtJ + J = DoF+... + DnO^F + O {O^+'F, d^) . (12) 

Here, we omitted all terms involving spatial derivatives. We introduced the 
transport coefficients Xm and Dm, 

x„ = $„(o) = (-0'" E ^. ' (0) ' 

™~ m! wi{0)---ujn{0) ■ ^ ' 
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If there is a clear separation of scales, i.e., Kn <C 1, it becomes possible 
to simplify the right-hand side of Eq. (12). In this case, DqF 3> DxdtF 
D2dfF, and higher derivatives of F can be neglected. However, in order 
to obtain a relaxation-type equation of motion for J, additional approxi- 
mations are required. For this purpose, one must impose the limiting pro- 
cedure, Xi — >■ 0, i > 2. In such a limiting procedure, all the poles of the 
retarded Green's function except the pole nearest to the origin (in the fol- 
lowing referred to as the "first pole") are pushed to infinity.^ Assuming that 
this can be done, we obtain the following equation of motion for J 

TRdtJ + J = DoF + DidtF + D2d^F + O {Dad^F, d^) , (15) 

where one can show that 



TR 



iwi(O)' 
Do = Gr{0,0) , 

= idu; Gr{u>,0) +Dotr. (16) 

ai=0 

Note that the time reversal symmetry of the retarded Green's function 
requires the first pole to lie on the imaginary axis, such that tr is real, 
as it should be. It is also important to remark that Eq. (15) describes the 
long- (but not asymptotically long-) time evolution of the system. This is 
the physical origin of the limiting procedure previously implemented and 
it is the only limit in which a relaxation equation can be obtained and, 
consequently, a relaxation time can be defined. If one wants to describe even 
shorter time scales, it becomes necessary to include more terms ~ Xi in the 
discussion. Then, the equations of motion will become more complicated. 
Naturally, the gradient expansion is recovered when all poles are pushed to 
infinity and, in this sense, the gradient expansion can be interpreted as an 
asymptotic solution of the more general equation of motion.^ 

Note that this limiting procedure is only allowed when the first pole lies 
on the imaginary axis. On the other hand, if the first pole and, for rea- 
sons of symmetry, its counterpart on the other side of the imaginary axis, 
have nonzero real parts one cannot disregard one of the poles while keeping 
the other. Then, the long-time evolution of the dissipative current will be 
characterized by a second-order differential equation and will display oscil- 
latory motion. Therefore, in this case the fluid-dynamical motion cannot 
be reduced to a simple relaxation equation, even in the small-frequency 
domain. 



January 14, 2013 17:38 



WSPC - Proceedings Trim Size: 9in x 6in Proceedings'Goa 



7 

6. Application: The Linectrized Boltzmann Equation 

The discussion presented above is general and can, in principle, be applied 
to any type of theory. In this section, we apply this formalism to kinetic 
theory. In particular, we calculate the shear viscosity and relaxation time 
coefficients for a weakly coupled classical gas of hard spheres via the Boltz- 
mann equation.'^ 

We perform all our calculations in the local rest frame of the fluid el- 
ement, where = (1,0,0,0). Since we restrict our discussion to shear 
viscosity, we consider the simplified scenario in which the fluid does not 
accelerate, u'^d^u" = 0, and does not expand, df^u^ = 0, and where tem- 
perature, T = f3Q^, and chemical potential, /x = ao/Po, are constant. With 
all these simplifications, the linearized Boltzmann equation can be written 
as^ 

dt6fu + yr-Vh-C6h=SiX,K), (17) 

where Sfu = fv - /ok, with /ok = exp (ao - /3o£^k), v = \<./E^, and S = 
/3o/oki^k'^^''fc''^'^M-Here, we introduced the collision operator, C, defined 
as 

X (^Sfp' + -^Sfp - -^^h' - -^^h) > (18) 

V/Op' /Op JOk' JOk / 

where M^kkz-s-pp/ is the transition rate, ly is the symmetry factor = 2 for 
identical particles) and dK = d'^k E^^ / {2tt)^ . 

Using Eqs. (17) and (18), it is possible to express the Fourier transform 
of the shear stress tensor, Tr^^'^{Q), in terms of the Fourier transform of the 
shear tensor, aa^{Q), in the form Tr*"^ = C^^^dap,^ where 

G^r^ (a;,q) = / dK k^^ k""^ ^ f3o fokE^^k<" k^> . (19) 

J iu! ~ iv ■ q — C 

In order to compute this Green's function, we define a function B°'^{Q, K) 
which satisfies 

{iuj-iwq-C) B''l^{Q,K)=S{Q,K)=(3oE^^k<"k^^fok. (20) 

As shown in the previous section, in order to extract the shear viscosity and 
relaxation time coefficients, it is sufficient to consider the case q = 0. Then, 
the dependence of B"^ on K can be expressed via the following expansion, 

oo 

B«'5(a;,i^) = /okfc<"fc'5>^a„(a;)£;£. (21) 

n=0 
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Substituting Eq. (21) into Eq. (19), we obtain the following expression for 
G^^"^(w,0), 



n=0 

oo 

= 2A''-«'3^7„+4,2a„(a;), (22) 



n=0 

where we introduced the thermodynamic integral 

= (2^TT)!! / ^^^oi^^r'' K - ElY . (23) 

Thus, the relation between -k^" and a'^" can be cast into a more con- 
venient form Tyi^'' (w, 0) = 2Gr (w, 0) a^^" (w, 0), with Gr (w, 0) being given 

by 

OO 

an(w). (24) 

n=0 

The whole frequency dependence of Gji{u),d) is contained in the func- 
tions a„(a;). The problem is then reduced to determining the analytic prop- 
erties of a„('^)- The way to solve this problem is to substitute the expansion 
(21) into Eq. (20), multiply by El^^k'^^' k , and integrate over dK. Then one 
obtains 

oo 

J2 (-iwl?"" + ^'"") an {OJ) = 0oIm+3,2 , (25) 

n=0 

where we defined the matrices 

^mn^f.ua0 = j dKE^^k^^" k^W {K) fokE]^k<" k>^\ (26) 

= / f^^^^^^ ("^' -Elf . (27) 
The solution of this equation can be formally expressed as, 

oo 

am(w) =/3o2Z [(-iwP + ^) 'J /„+3,2. (28) 

Tl=0 

Thus, we obtain 



-I mn 



U-iujV + A) M /„+3,2. (29) 



m=0 n=0 
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The poles of this function can be obtained from the roots of the determinant 
det {—i(jjV + A) =0. Thus, the relaxation times and viscosity coefficients 
are determined by inverting and finding eigenvalues of infinite-dimensional 
matrices. In practice, however, the expansion (21) has to be truncated and 
one never actually deals with infinite matrices. This is a common approach 
and happens quite often when dealing with the collision operator. Let us 
try the simplest case and consider only one term in the expansion (this 
corresponds to using the Israel- Stewart 14-moment approximation in the 
moments method) . Then all matrices become numbers and everything sim- 
plifies considerably. In this case, the retarded Green's function has only one 

polo, Wo = -M°7/42.^ 

Since in this case the thermodynamic force is given by _F = 2a^^ , the 
shear viscosity coefficient is given by Dq. Thus, the shear viscosity and 
relaxation time are given by 



This implies that ry/rvr = /3oJ32. In the masslcss limit, for a gas of hard 
spheres, = (3/5)142 nocr, with a being the total cross-section. 

Exactly the same results were obtained in Ref.^^ within a completely 
different approach. This clearly demonstrates that the relaxation time in 
Israel- Stewart theories, which determines the time scale related to transient 
dynamics, is indeed a microscopic scale. In other words, the relaxation time 
is determined by the inter-particle scattering rate and not by an arbitrary 
fluid-dynamical time scale. 

7. Conclusion 

In this work we proved that the existence of singularities in the retarded 
Green's functions considerably restricts the applicability of the gradient ex- 
pansion as a method to extract macroscopic dynamics. As a matter of fact, 
the causal structure of relativistic fluid dynamics hinges on the presence of 
poles in the Green's functions. 

We computed the first pole of the retarded Green's function associated 
with shear stress in a weakly coupled classical gas of hard spheres. We 
proved that the inverse of this pole is directly related to the relaxation 
time that appears in Israel-Stewart theory. Consequently, the relaxation 
time in Israel-Stewart theories corresponds to a microscopic time scale and 
such theories are indeed describing the real transient dynamics of the fluid. 




I42I32 



(30) 



(31) 
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In summary, the true relaxation time is always given by the first pole of 
the retarded Green's function. In general, the location of this pole cannot 
be found using a truncated Taylor expansion around the origin or, in other 
words, via the gradient expansion. This result should be expected since the 
transient behavior of a given system cannot be unambiguously extracted 
from its asymptotic solution given by the gradient expansion. 
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insightful discussions. H.N. was supported by the ExtrcMe Matter Institute 
EMMI. The authors thank the Helmholtz International Center for FAIR 
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